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We investigate two complementary problems related to maintaining the relative positions of N 
random walks on the line: (i) the leader problem, that is, the probability £at(£) that the leftmost 
particle remains the leftmost as a function of time and (ii) the laggard problem, the probability 
IZNit) that the rightmost particle never becomes the leftmost. We map these ordering problems 
onto an equivalent (N — l)-dimensional electrostatic problem. From this construction we obtain a 
very accurate estimate for £jv(t) for N = 4, the first case that is not exactly soluble: C^it) oc £ _/34 , 
with p4 = 0.91342(8). The probability of being the laggard also decays algebraically, IZNit) oc t~ lN ; 
we derive 72 = 1/2,73 = 3/8, and argue that 7jv — > iV -1 In N as N — > 00 . 

PACS numbers: 02.50.Cw, 05.40.-a, 05.50.+q, 87.18.Sn 



I. INTRODUCTION 

Consider N identical and independent random walkers 
that are initially located at xi(0) < £2(0) < • * * < %n(Q) 
on an infinite line Q. There are many interesting ques- 
tions that can be posed about the order of the particles. 
For example, what is the probability that all walkers 
maintain their relative positions up to time t, that is, 
X\(t) < X2(t) < • • • < xn(t) for all < r < tl By map- 
ping this vicious random walk problem onto the diffusion 
of a single effective particle in TV dimensions and then ex- 
ploiting the image method for the diffusion equation, this 
ordering probability was found |) to decay asymptot- 
ically as t~ aN with a N = N(N - l)/4. Many additional 
aspects of this problem have been investigated within the 
rubrics of vicious random walks @, §, f?|, H |§] and 
"friendly" random walks |l(], |ll| . 

In this work, we study two related and complemen- 
tary random walk ordering problems. In the "leader" 
problem, we ask for the probability C^{t) that the ini- 
tially leftmost particle in a group of N particles re- 
mains to the left of all the other particles up to time 
t 0, [OJ [H| 0. In the "laggard" problem, 
we are concerned with the probability IZNif) that the 
initially rightmost particle from a group of N particles 
never attains the lead (becomes leftmost). These two 
probabilities Cn(P) and 7Z]y(t) decay algebraically 



Cjsf(t) oc t 



7£/v(t) (X t 



(i) 



as t — > 00 and our basic goal is to determine the expo- 
nents Pn and 7at. 

These ordering problems arise in a variety of contexts. 
Physical applications include wetting phenomena |2], [j| 
and three-dimensional Lorentzian gravity fig]. A more 
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probabilistic application is the ballot problem |T^] , where 
one is interested in the probability that the vote for a 
single candidate remains ahead of all the other candidates 
throughout the counting; this is just a restatement of the 
leader problem. Another example is that of the lamb and 
the lions [l5|, |l6| , in which one is interested in the survival 
of a diffusing lamb in the presence of many diffusing lions. 
In one dimension, a lamb that was initially in the lead 
must remain the leader to survive. 

For the leader problem, exact results are known for 
small N only: fa = 1/2 and fa = 3/ 4 i, O, |o), 
while (3 N In (AN)/ 4 for large TV (l|, gj, |16|, |l7j. This 
slow increase arises because adding another particle has 
little effect on the survival of the leader when N is large. 
For N > 4, no exact results are available and one focus 
of our work is to obtain an accurate estimate of (3n for 
the case N = 4. We accomplish this by mapping the re- 
action onto an equivalent electrostatic potential problem 
due to a point charge within an appropriately-defined 
three-dimensional domain. This mapping provides both 
an appealing way to visualize the reaction process and 
an accurate estimate of the survival exponent fa. 

For the laggard problem, we employ the same method 
as in the leader problem to obtain 73 = 3/8. We 
also estimate the asymptotic behavior of 1Z n and find 
7at — > N~ 1 lnN as N — > 00. As is expected, a laggard 
in a large population likely remains a laggard. Therefore 
the probability of remaining a laggard decays very slowly 
with time for large N. 

In the next section, we review known results about the 
leader problem in a 3-particle system. In Sec. Ill, we 
outline the electrostatic formulation of the leader prob- 
lem and then apply it to the 4-particle system in Sec. IV. 
A numerical solution of the pertinent Laplace equation 
gives fa = 0.91342(8), a significant improvement over 
the previously-quoted estimate fa w 0.91 Q. In Sec. V, 
we turn to the laggard problem and give an asymptotic 
estimate for the exponent 7at. Concluding remarks and 
some open questions are given in Sec. VI. 
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II. CONVENTIONAL APPROACH TO THE 
LEADER PROBLEM 

We begin by mapping the original problem of N dif- 
fusing particles xi, x<i, . . . , xn on the line onto a single 
effective diffusing particle located at (xi, X2, . . . , xjy) in 
TV-dimensional space. The particle order constraints on 
the line translate to bounding hyperplanes within which 
the effective particle is confined [[| [2(J . The effective par- 
ticle is absorbed if it hits one of these boundaries. For 
the leader problem, the explicit shape of these bound- 
ing hyperplanes can be easily worked out for the cases 
N = 2 and N = 3; we will later extend this analysis to 
the 4-particle system. 

For 2 particles, their separation undergoes simple dif- 
fusion and the process terminates when the separation 
equals zero. Thus the survival probability of the leader 
decays as t -1 / 2 . To fix notation and ideas for later sec- 
tions, we now study the 3-particle system. For a leader 
at x\{t) and particles at ^(t) and xs(t), we view these 
coordinates as equivalent to the isotropic diffusion of a 
single effective particle at (xi(i), X2(t), x% (£)) in three di- 
mensions. Whenever this effective particle crosses the 
plane : X{ = xj, the original walkers at xi and Xj 
in one dimension have reversed their order. There are 
three such planes Ai2,Ai3,A23 that divide space into 6 
domains, corresponding to the 3! possible orderings of 
the three walkers (Fig. |l|(a)). These planes all intersect 
along the (1,1,1) axis. 

We may simplify this description by projecting onto 
the plane x\ + X2 + xs =0 that contains the origin and 
is perpendicular to the (1,1,1) axis. Now the plane A±2 
may be written parametrically as (a, a, b) and its intersec- 
tion with the plane X1+X2+X3 = is the line (a, a, —2a). 
Likewise, the intersections of A13 and A23 with the plane 
x\ + X2 + X3 = are (a, —2a, a) and (—2a, a, a), respec- 
tively (Figs. 0(b)). 

The survival of the leader corresponds to the effective 
particle remaining within the adjacent domains 123 and 
132 in Fig. 0(b). The background particles at X2 and X3 
are allowed to cross, but the leader at x\ always remains 
to the left of both X2 and X3. The union of these two 
domains defines a wedge of opening angle 120°. Since 
the survival probability of a diffusing particle within a 
wedge of arbitrary opening angle (p and absorbing bound- 
aries decays asymptotically as t _7r / 2< ^ pl| , we deduce the 
known result that the leader survival probability expo- 
nent is /?3 = 3/4. 



III. ELECTROSTATIC FORMULATION 

For more than 3 particles, the domain for the effec- 
tive particle is geometrically more complex and the cor- 
responding solution to the diffusion equation does not 
seem tractable. We therefore recast the survival prob- 
ability of the effective diffusing particle in terms of the 
simpler problem of the electrostatic potential of a point 




FIG. 1: The order domains for 3 particles in (a) the full 3- 
dimensional space, and (b) projected onto the subspace per- 
pendicular to the (1,1,1) axis. The notation ijk is short- 
hand for Xi < Xj < Xk- The allowed region corresponding to 
survival of a leader (x± < ^25X3) is indicated by the darker 
shading, while the lighter shaded region corresponds to the 
laggard problem (#3 it xi, re- 
charge within the same geometric domain [^0|. Let S(t) 
be the survival probability of a diffusing particle within 
an infinite wedge-shaped d-dimensional domain with ab- 
sorbing boundary conditions. Let V(r) be the electro- 
static potential due to a point charge within the same 
domain, with V = on the boundary. Generically, these 
two quantities have the asymptotic behaviors 

S(t) oc r$ , t -> 00 , 

V(r) ex r _/i , r -> 00 . (2) 

More relevant for our purposes, these two quantities 
are simply related by Jl^, [l^, |20| 

rt fVt 
J S(t)dt~ J V^r^dr. 

This equivalence arises because the integral of the dif- 
fusion equation over all time is just the Laplace equa- 
tion. Thus the time integral of the survival probability 
has the same asymptotic behavior as the spatial integral 
of the electrostatic potential over the portion of the do- 
main that is accessible by a diffusing particle up to time 
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t. Substituting the respective asymptotic behaviors from 
Eqs. (H), and noting that the allowed wedge domain for 
an TV-particle system has dimension N — 1, we obtain the 
fundamental exponent relation 



Thus the large-distance behavior of the electrostatic po- 
tential within a specified domain with Dirichlet boundary 
conditions also gives the long-time survival probability 
of a diffusing particle within this domain subject to the 
same absorbing boundary conditions. From this survival 
probability, we can then determine the original ordering 
probability. 

To illustrate this approach, let us determine the various 
ordering probabilities of 3 particles on the line in terms 
of the equivalent electrostatic problem. In fact, it is sim- 
pler to work backwards and find the equivalent ordering 
problem that corresponds to a specific wedge domain, For 
example, consider the 60° wedge 123 in Fig. 0. If the ef- 
fective particle remains within this wedge, the initial par- 
ticle ordering on the line is preserved. This corresponds 
to the vicious random walk problem in which no particle 
crossings are allowed. To obtain the asymptotic behavior 
of the potential of a point charge interior to this wedge, 
let us assume that V(r) ~ r _/x f((f) as r — > oo. Substitut- 
ing this ansatz into the 2-dimensional Laplace equation, 
we obtain the eigenvalue equation f"(cp) = — /i 2 /((^), 
subject to f((f>) = on the wedge boundaries. For the 
60° wedge, the solution with the smallest eigenvalue is 
f((p) = sin(3(/?). Thus /i = 3, leading to the known result 
(3 = (fi — 7V + 3)/2 = 3/2, for 3 vicious walkers. 

IV. 4-PARTICLE SYSTEM 

The state of the system may be represented by an ef- 
fective diffusing particle in 4 dimensions. By projecting 
onto the 3-dimensional subspace x\ + X2 + £3 + £4 = 
that is orthogonal to the (1, 1, 1, 1) axis, the order do- 
mains of the original particles can be reduced to 3 spatial 
dimensions. In this 3-dimensional subspace, the bound- 
ary A12 : x\ = X2 becomes the plane (a, a, 6, c), with 
2a + b + c = 0. Likewise, A13 may be written paramet- 
rically as (a, 6, a, c), where again 2a + b + c — 0. Thus 
the locus L123 = A\2 H A13 : x\ = X2 = £3 is the line 
(a, a, a, —3a). This body diagonal joins the nodes 4 and 
4 in Fig. ||. Along this axis, the original particle coor- 
dinates on the line obey the constraint X\ = X2 = £3, 
with X4 < xs on the half-axis closer to node 4 and 
X4 > xs on the half-axis closer to node 4. A similar 
description applies to the axes L124 = (a, a, —3a, a) be- 
tween 3 and 3, L134 = (a, —3a, a, a) between 2 and 2, and 
£234 = (—3a, a, a, a) between 1 and 1. 

The locus where x\ = X2 and x% = £4 simultaneously, 
is the line £12,34 = A12 D A34 = (a^a^ —a^ —a). Like- 
wise, £13,24 = — fl, fl, — and £14,23 = —a, —a, a). 
Viewed in the orthogonal 3-subspace, the 6 planes A^ 




14,23 



FIG. 2: Top: Cardboard model of the ordering domains for 
4 particles on the line after projection into the subspace per- 
pendicular to (1,1,1,1). This structure consists of only 6 
intersecting planes. Each plane bisects the cube and is de- 
fined by the equality of two coordinates. Lower: Schematic 
of the same system. The wedge ijkl denotes the region where 
Xi < Xj < Xk < xi. The union of six such ordering wedges, 
corresponding to the leader problem x\ < X2,X3,X4, are la- 
beled. This domain is bounded by the rays 02, 03, 04, where 
the origin is the dashed circle at the intersection of the axes 
£14,23 = u, £12,34 = v, and £13,24 = w. One constraint plane, 
x\ = £3, is shown shaded (outside the cube). 



intersect in 7 lines (4 and 3 Lij^i), and divide 

the subspace into 24 semi-infinite wedges, as shown in 
Fig. cube(b). Each of these wedges corresponds to one 
of the 4! orderings of the walkers in one dimension. 

We first illustrate the electrostatic formulation of this 
system for the vicious random walk problem. Since 
the initial particle order is preserved, the effective dif- 
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fusing particle remains within a single wedge ijkl in 
Fig. g. As outlined in the previous section, the sur- 
vival probability of this effective particle corresponds to 
the electrostatic potential of a point charge within this 
one wedge, with the boundary surfaces held at zero po- 
tential. To solve this electrostatic problem, it is conve- 
nient to place the point charge at the symmetric loca- 
tion (u,v,w) = (0, 1, 1/2) within the wedge 1234, where 
the (u,v,w) axes are defined in Fig. |2| From the image 
method, the potential due to a point charge within one 
wedge is equivalent to the potential of an array of 24 
symmetrically-placed point charges consisting of the ini- 
tial charge and 23 image charges, with positive images at 
-(0,1,1/2), ±(0,1,-1/2), ±(±1/2,0,1), ±(1, ±1/2,0), 
and negative images at ±(±1/2, 1, 0), ±(0, ±1/2, 1), and 
±(0,1, ±1/2). Using Mathematica the asymptotic be- 
havior of the potential in wedge 1234 (where the original 
charge is placed) due to this charge array is 



V(r) 



a 2 r 



-li 



a 3 r 



-13 



-15 



(The coefficients ai depend on the location of the charge 
and on the orientation of r.) Using the exponent rela- 
tion (||), the asymptotic survival probability of 4 vicious 
walkers is given by 



S(t) = M~ 3 ± b 2 t~ 5 ± b s r 



As 



The leading behavior confirms the known result 
and we obtain the form of the corrections as wel 
an amusing aside, notice that an r~ 7 dependence for the 
potential is normally achieved by a 64-pole charge con- 
figuration. The high symmetry of the 24 charges in the 
ordering problem leads to a multipole field normally as- 
sociated with at least 64 point charges. 

We next turn to the leader problem. This system cor- 
responds to the electrostatic problem within the com- 
bined domain of the 6 wedges marked 1234, 1243, 1423, 
1432, 1342, and 1324, in Fig. |. The resulting domain is 
a tetrahedral corner, with its apex at the center of the 
cube, that is flanked by the rays 02, 03, and 04. Despite 
the simplicity of this domain, we are unable to solve this 
electrostatic problem analytically and we have instead 
studied the problem numerically. 

In the allowed region of the cube in Fig. |[ we dis- 
cretize space and solve the electrostatic potential of a 
point charge by using successive over-relaxation, with 
the domain boundaries at zero potential. For simplic- 
ity, the charge is chosen to be at the symmetric point 
(1,1,1). For the outer faces of the cube we use two dif- 
ferent boundary conditions: (a) absorbing (V = 0), and 
(b) reflecting, (dV/dn = 0, where n is the direction nor- 
mal to the surface). The true potential - that of the in- 
finite wedge - lies between these two extremes. We also 
exploit the symmetry about the (1,1, 1) diagonal and use 
only the domain marked 1234 in Fig. g, with absorbing 
boundary conditions on the plane 3434, and reflecting 
boundary conditions on the planes 1212 and 1414. As al- 
ready discussed, the outer cube face, 1324, is taken to be 



reflecting or absorbing. This space savings allows us to 
carry out computations for a cube of 500 lattice spacings 
on a side. 



2.827 
2.826 
2.825 
2.824 
2.823 
2.822 
2.821 
2.820 



V 




0.02 



0.04 0.06 
1/r 



0.08 



0.1 



FIG. 3: Local exponent /i(r), as a function of 1/r for the 
tetrahedral wedge with absorbing (+) and reflecting (o) finite- 
size boundary conditions at the outer face of the cube. 

While the power-law decay of the potential sets in 
quickly, the finite-size effect is pronounced and it is per- 
ceptible already at 25 lattice spacings away from the 
charge. This is the primary limitation on the accu- 
racy of our exponent estimate. Nevertheless, the lo- 
cal exponent ji(r) = —d\nV(r)/d\nr varies only at the 
fourth digit (Fig. ^). The approach of the local expo- 
nent to the asymptotic limit also suggests that V(r) 
has the form V(r) ~ r _/i ± Ar~ 4 . Assuming that this 
is the case, extrapolation of the data in Fig. || gives 
\i = 2.82684 ± 0.00016, where the error bar is the dif- 
ference in the extrapolated value of ja(r) from the two 
different boundary conditions. From the exponent cor- 
respondence given in Eq. (||), we thereby obtain, for the 
lead probability, 



U(t) ~ ± At~ 3/2 , /? 4 = 0.91342(8) 



(4) 



It is hard to match this numerical accuracy with that 
from direct simulations of the survival of the leader. We 
simulated 10 9 realizations of the system; this gives ex- 
tremely linear data for the time dependence of the leader 
survival probability on a double logarithmic scale. To 
estimate the exponent /?4, we computed the local slopes 
of the survival probability versus time in contiguous time 
ranges between t and l.bt when plotted on a double log- 
arithmic scale. These local exponents are plotted against 
l/lni.st (Fig. |^). The results are compatible but much 
less accurate than Eq. (0). 



V. THE LAGGARD PROBLEM 

In the laggard problem, we study the probability that 
the initially rightmost particle at xn has never been the 
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FIG. 4: Direct simulation results for the local exponent in 
the survival probability for 10 9 configurations. The arrow 
indicates our estimate for (3a from Eq. (W). 



leader during the time interval (0,£). The laggard prob- 
lem can also be recast into the diffusion of a single ef- 
fective particle within an TV-dimensional wedge-shaped 
region, with absorbing domain boundaries. This map- 
ping leads to the basic conclusion that every particle that 
is initially not in the lead exhibits the same asymptotic 
behavior as the last particle. Indeed, for any particle i 
initially at Xi, the regions X{ x±, . . . , Xi-i, Xi+i, . . . , xn 
are isomorphic. The initial condition merely fixes the lo- 
cation of the effective particle in this allowed region. An- 
other fundamental observation is that the allowed regions 
of the effective particle for the leader and the laggard 
problems are complementary for all N. 

For two particles, the probability that the laggard does 



-1/2 



not become the leader obviously decays as t 
72 = 1/2. The case TV = 3 is more interesting but also 
solvable. The condition that a particle initially at xs 
never attaining the lead (xs £1,^2) is equivalent to 
the effective particle remaining with the lighter shaded 
region in Fig. 0(b). Since the opening angle of the re- 
sultant wedge is ip = 47r/3, the corresponding survival 
probability decays as t _3//8 , implying that 73 = 3/8. 

We have performed direct numerical simulations of the 
process to estimate the exponent 77V for N = 3,4,5, 
and 6. Each simulation is based on 10 6 realizations in 
which N — 1 particles are initially at the origin, while 
the laggard is at x = 1. Each simulation is run until the 
laggard achieves the lead or 10 5 time steps, whichever 
comes first. From the survival probability, we estimate 
7at ~ 0.35,0.30,0.26, and 0.23 for N = 3 - 6, respec- 
tively. Since we know that 73 = 0.375, the discrepancy 
of 0.025 between the simulation result and the theory is 
indicative of the magnitude of systematic errors in this 
straightforward numerical approach. 

While it appears difficult to determine the exponent 7 at 
analytically for general N > 3, the situation simplifies 



in the large TV limit because the position of the leader 
becomes progressively more deterministic. Indeed, the 
probability density Pjy(y,t) that the leader is located at 
distance y from the origin (assuming that all particles are 
initially at the origin) is Jl6[ ] 

Ne -y 2 /4Dt r 1 / M^- 1 




where erfc(x) = 1 — erf(x) is the complementary error 
function. Performing a large N analysis, we find that 
the probability density Pjsr(y,t) approaches a Gaussian 



PN(v,t) 



1 



exp 



f (y + y*) 2 \ 



(5) 



with the mean and the variance of this distribution given 
by 



y* 



Dt 
72" 



(6) 



where z is determined from the transcendental relation 

N 



ze = 



4tt 



(7) 



Consequently, the parameter z diverges as z « Vln N 
when N — > 00. The ratio of the dispersion to the mean 
displacement thus decreases as cr/y* « (21niV) _1 , so 
that the position of the leader indeed becomes more de- 
terministic as TV — » 00. 

Therefore in the large N limit we can assume that the 
leader is moving deterministically and its position is given 
by —y*(t). Then the probability that the laggard never 
achieves the lead is equivalent to the probability that a 
diffusing particle initially at the origin does not overtake 
a receding particle whose position is varying as —y*(t). 
This corresponds to the solution to the diffusion equa- 
tion in the expanding region x G (—?/*(£), 00) with an 
absorbing boundary condition at the receding boundary 
x = —y*. When y* ex t 1 / 2 , this diffusion equation can 
be solved exactly by reducing it to a parabolic cylinder 
equation in a fixed region. However, in the limit t —> 00, 
we can obtain asymptotically correct results much more 
simply. Because the absorbing boundary recedes from 
the laggard particle relatively quickly, we solve the prob- 
lem by assuming that the density P{x,t) of the laggard 
particle approaches a Gaussian with yet unknown ampli- 
tude ll N (t) 



P(x,t) 



n N (t) 



exp 



x 2 \ 
~4Dt J ' 



(8) 



Although this distribution does not satisfy the absorbing 
boundary condition, the inconsistency is negligible since 
the exponential term in Eq. (||) is of order iV _1 at the 
boundary y = y*. 

The probability 1ZN(t) is now found self-consistently 
by equating the "mass" loss to the flux: 



oIIZn 
dt 



= D 



dP 

dx 



(9) 
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Using Eq. (m to compute the flux we convert Eq. (ph into 



dll 



N 



dt 



z 1 
N 



(10) 



from which the exponent 7at is z 2 /N. This is of course 
valid only in the large N limit. Taking this limit in 
Eq. (]?]) we obtain 



IniV lnlnTV 



7iv 



N 



27V 



(11) 



Thus as N gets large, the probability that the laggard 
never attains the lead decays extremely slowly with time. 
This fits with the naive intuition that if the number of 
particles is large a laggard initially is very likely to remain 
a laggard. Each additional particle makes it even less 
likely that the laggard could achieve the lead. Amusingly, 
this asymptotic exponent predictions is numerically close 
to the previously-quoted results from direct numerical 
simulations of the laggard problem. 



VI. CONCLUDING REMARKS 

We investigated two dual random walk ordering prob- 
lems in one dimension: (i) what is the probability that a 
particle, that is initially in the lead, remains in the lead 
and (ii) what is the probability that a particle, that is 
initially not in the lead, never achieves the lead? These 
problems are most interesting in one spatial dimension 
because of the effective correlations between the inter- 
acting particles. These correlations are absent in two 
dimensions and greater, so that an TV-particle system re- 
duces to a 2-particle system |l5| (Hf . 

We determined the respective exponents (3n and 7at 
associated with the lead and laggard probabilities for 
general TV. Both exponents can be determined by ele- 
mentary geometric methods for N = 2 and 3 and by 



asymptotic arguments for N — > oo. Our new results are 
the following: (i) a precise estimate for (3^ and (ii) the 
large- N behavior of 7at. 

A simple generalization is to allow each particle i to 
have a distinct diffusion coefficient Di. The exponents 
Pn and 7at will now depend on the diffusion coefficients, 
except for N = 2, where fa and 72 always equal 1/2. 
The case N — 3 is still solvable by introducing rescaled 
coordinates yi = Xi/y/Dl to render the diffusion of the 
effective particle isotropic, after which the mapping to 
the wedge can be performed straightforwardly. We thus 
find 



ft 



73 



2 cos 

7T 



COS 



^(D 1 ^D 2 )(D 1 ^D 3 ) 



An amusing special case is the case of a stationary lag- 
gard, for which 73 = 1/3. 

Finally, it is also worth mentioning a promising de- 
velopment to solve the diffusion equation in the domains 
defined by the ordering of one-dimensional random walks. 
This is the recent discovery of deep connections between 
vicious walkers and random matrix theory j(| ^, ^ . These 
allow one to not only re-derive the exponent a n of the 
original vicious random walk problem, but also lead to 
many new results. It would be extremely useful if these 
techniques could be extended to the leader and laggard 
problems. 
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